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Abstract. The factorization of the amplitude of the process 7*7 — » 77 in the low energy 
and high photon virtuality region is demonstrated at the Born order and in the leading loga- 
rithmic approximation. The leading order two photon (generalized) distribution amplitudes 
exhibit a characteristic InQ 2 behaviour and obey new inhomogeneous evolution equations. 
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1. INTRODUCTION 

The pointlike coupling to quarks of the photon enables to calculate perturbatively 
part of its wave function. A twist expansion generates non-leading components of 
the photon distribution amplitude [TJ, from which the lowest order one is chiral- 
odd and proportionnal to the magnetic susceptibility of the vacuum. The study of 
the two photon state is kinematically richer and is thus a most welcome theoretical 
laboratory for the study of exclusive hard reactions. 

The parton content of the photon has been the subject of many studies since the 
seminal paper by Witten [2J. A recent paper [3] extended the notion of anomalous 
parton distribution in a photon to the case of generalized parton distributions 
(GPDs) used for the factorized description of the non-diagonal kinematics of deeply 
virtual Compton scattering (DVCS) on a real photon target, 7* (9)7 — > 77, namely 
at large energy and small hadronic momentum transfer but large photon virtuality 
(Q 2 = -q 2 )- 

As for the two meson case jl] the two photon generalized distribution amplitudes 
describe the coupling of a quark antiquark (or gluon-gluon) pair to a pair of 
photons, and are related by crossing to the photon GPDs. 

We study [6] the scattering amplitude of the 7* (5)7 — > 77 process in the near 
threshold kinematics, namely at small s and large — t ~ Q 2 , at large Q 2 and in the 
leading order of the electromagnetic and strong couplings. This enables us to define 
and calculate perturbatively the Born approximation of the diphoton GDAs. 

This is reminiscent of, but quite different from, the perturbative calculation of 
the rho rho GDA in terms of the rho DA [7] , for which both incoming photon were 
chosen to be hard in order to justify this factorization of the GDA. 



2. 7*fah(Pi) 7(^)7(^2) NEAR THRESHOLD 



Two photon production in Compton scattering on a photon target 

7*(?)7(?0 ->7(Pi)7(P2) (1) 

involves, at leading order in a em , and zeroth order in 0:5 six "box" Feynman 
diagrams with quarks in the loop. 

Restricting to the threshold kinematics where W 2 = (P1+P2) 2 = simplifies 
greatly the tensorial structure of the amplitude while still preserving the richness 
of the skewedness (£) dependence of GDAs. Our conventions for the kinematics are 
the following: 

Q 2 , Q 2 f0 . 

q = p n , q = — n , (2) 

s s 



Pi = (P , P2 = (p , C = 1 - C» ( 3 ) 

where p and n are two light-cone vectors and 2p ■ n = s. The momentum I in the 
quark loop is parametrized as 

l» = ZI f + {3 n » + l T , (4) 
with l\ = — l 2 . The scattering amplitude is written as 

^ = e^e^T^, (5) 

where the four photon polarization vectors are transverse with respect to Sudakov 
vectors p and n. 

The tensorial decomposition of T MiyQ/3 (Vy = 0) reads 

= \g^gfW x + I + g^gf - g^gf) W 2 + \ - <tffi>) W 3 , 

(6) 

and it involves three scalar functions Wi, i = 1,2,3. 

The integration over / is performed as usual within the Sudakov representation, 
using 

d 4 l = -dzd(3d 2 lT^ — dz df3 d\ 2 . 

Let us note that in order to interpret our result in terms of factorized quantities, 
we will keep our expressions unintegrated with respect to z, the mass of the quark 
playing the role of an infrared regulator. 

One first integrates in j3 using the Cauchy theorem. The propagators induce 
poles in the complex /5-plane and the pole positions depend on the values of z and 
£. The four poles lie all below the real axis for z > 1 and lie all above the real axis 
for z < 0, the only region where the amplitude may not vanish is 1 > z > 0. This 
leads to a natural interpretation of partonic fraction of momentum. One 

then identifies different regions defined from the relative values of z, ( and 1 — (. 



This is reminiscent of the different regions encountered in the kinematics of the 
generalized parton distributions H(x,£,t), with the boundaries controlled by the 
relative values of x and £, where x and £ are the quantities related to our z and ( 
variables. For each of the six contributing diagrams the remaining integral over l 2 
contains a UV divergent part which cancels in their sum. We get 
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Let us now interpret the results ^ and (|9| from the point of view of QCD 
factorization based on the operator product expansion. 



3. QCD FACTORIZATION AND THE 77 GDA 



Let us consider two quark non local correlators on the light cone and their matrix 
elements between the vacuum and a diphoton state which define the diphoton GDA 



F q 



d V A 



1(22-1)1 



2tt 



{l{Pi)l{P2)\q{ 



-yN 



-gTe;( Pl )e:(p 2 )^(z,C,0) 

(10) 

where we denote N = n/n.p and where we did not write explicitely the electromag- 
netic and the gluonic Wilson lines. We need also to define the matrix element of 
the photonic correlator 



dy 



27r ^- 1 )i( T (p 1 ) T (p 2 )| J P^(-|iV)if(|iV)|0) 



(11) 



where F N ^ = N V F U » , which mixes with the quark correlator rtlOh although they 
are not of the same order in a em [2J. 

We regulate through the usual dimensional regularization procedure the UV 
divergent quark correlator matrix elements and obtain (with 4 = - +7£ — log47r) 



N c e 2 



+ logm 2 



(12) 



FIGURE 1. The unpolarized anomalous diphoton GDA §\/(N c e 2 q /{2ir 2 ) log ^) at 
Born order and at threshold for £ = 0.1 (dashed), 0.2 (dash-dotted), 0.4 (solid). 



with F(z,() = 

M J( ,_ + + £Pi^W) , K _ + fP^lzO , (f _ . 

(13) 

The ultraviolet divergent parts are removed through the renormalization proce- 
dure (see for example [8]) involving quark and photon correlators (O q , O 7 ) corre- 
sponding to q(-$N)'y.Nq(%N) and F NfJ, (-^N)F^(^N). The renormalized oper- 
ators are defined as : 

(14) 

The matrix element of the renormalized quark-quark correlator is thus equal to 

< 7(Pi)7(P2)|O&|0 >= Z qq < i{pMp2)\O q \0 > +Z qi < 7(pi)7(p 2 )|O 7 |0 > (15) 
with Z qq = 1 + fyj . Since the matrix element < 7(^1)7(^2) \O q \0 > contains a UV 




divergence (Eqn. 12) and since < 7(^1)7(^2) |O 7 |0 > is UV finite and of order a° m , 
one can absorb this divergence into the renormalization constant Z qi . The normal- 
ization of the renormalized correlator is fixed with the help of the renormalization 
condition which is chosen as < 7(^1)7(^2) |C^.|0 >= at the renormalization scale 
Mr = m. In this way the renormalized GDA is equal to 

N r e 2 m 2 
K. = -^g^e^ Pi y v {p 2 )\og^F{zX) • (16) 

R 

from which we obtain - after identifying the renormalized scale with the factoriza- 
tion scale, M R = M F - that $l(z,(,0) = -^logf£-F(z,()- This expression to- 
gether with the Born order coefficient function C v = e 2 Q — |) leads to the quark 
contribution to the 7*7 — > 77 scattering amplitude 

1 

W! = J dzC v (z)$l(z,(,0) ■ (17) 



FIGURE 2. The polarized anomalous diphoton GDA $l/{Nc e 2 /(27r 2 ) log ^L) at Born 
order and for ( = 0.1 (dashed), 0.2 (dash-dotted), 0.4 (solid). 



The contribution to W\ in Eq. ([7]) related to the photon operator involves a 
new coefficient function of order a 2 em calculated at the factorization scale Mp 
and convoluted with the zeroth order in a em part of the photon GDA. This 
contribution coincides with the expression ^ in which the quark mass m is 
replaced by the factorization scale Mp playing now the role of infra-red cutoff. The 



sum of contributions related to the photon operator and to the quark operator (17) 
reproduces then the result of direct calculations (|7h. Let us note, that by choosing 
the factorization scale Mp equal to the hard scale of our process Q, Mp = Q, the 
resulting W\ comes only from the quark GDAcontribution. 

We have thus demonstrated that it is legitimate to define the Born order 
diphoton GDAs at zero W as 
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The same procedure can be applied to the axial vector correlator defining $ 3 

d V i(2z-m /../_ \ _ ./„ \ i — / y»A_. at.. _,V 



We get 
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Since we focused on the logarithmic factors, we only obtained the anomalous part 
of these GDAs. Their z— and £ — dependence are shown on Figs. [T] and |2j Note that 
they are discontinuous functions of z at the points z = ( and z = (. 
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